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Abstract. This article presents an introduction to optical pumping, 
atomic polarization and the Hanle effect in weakly magnetized stellar 
atmospheres. Although the physical processes and the theoretical frame- 
work described here are of interest for applications in a variety of astro- 
physical contexts (e.g. scattering polarization in circumstellar envelopes 
and polarization in astronomical masers), the article focuses mainly on 
the quest for understanding the physical origin of the linearly polarized 
solar limb spectrum. It considers also the development of the Hanle ef- 
fect as a reliable diagnostic tool for making feasible new advances in solar 
photospheric and chromospheric magnetism. Particular emphasis is given 
to a rigorous modeling of polarization phenomena as the essential link be- 
tween theory and observations. Some of the most recent advances in this 
field are presented after carefully explaining how the various radiation 
pumping mechanisms lead to atomic polarization in the absence and in 
the presence of weak magnetic fields. 



1. Introduction 

Probably the first thing I should point out is that the "second solar spectrum" is 
nothing but the observational signature of the "order" that exists in the atomic 
system^. This "atomic organization" is what we call atomic polarization (i.e. 
the existence of population imbalances among the sublevels of any given degen- 
erate atomic or molecular level and/or the presence of quantum interferences or 
coherences between any given pair of sublevels, even among those pertaining to 
different levels). But, what is forcing the ions, atoms and molecules of the stel- 



1 Consejo Superior de Investigaciones Cientfficas, Spain. 

2 The "second solar spectrum" is a term adopted by Stenflo and Keller (1997) to refer to the 
linearly polarized solar limb spectrum which can be observed with spectropolarimeters that 
allow the detection of very low amplitude polarization signals (with Q/I of the order of 10 -3 
or smaller) . Weak polarization signals in spectral lines had previously been observed in promi- 
nences outside the solar limb (e.g. Lyot 1934; Hyder 1965; see also the review by Leroy 1989) 
and on the solar disc close to the limb (e.g. Redman 1941; Bruckner 1963; Wiehr 1978; Stenflo 
et al. 1983), but most of the structural richness of the linearly polarized spectrum had re- 
mained inaccessible. The observations of Stenflo and coworkers with the polarimeter ZIMPOL 
(see also the atlas of Gandorfer, 2000) have been confirmed (and extended to the full Stokes 
vector) by Dittmann et al. (2001), Martinez Pillet et al. (2001) and Trujillo Bueno et al. 
(2001) using the Canary Islands telescopes. 
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lar atmospheric plasma to behave this way? As we shall see below, the atomic 
polarization is the result of a transfer process of "order" from the radiation field 
to the atomic system. It is thus natural that the logical structure of this review 
article is as follows: quantification of the "order" of the radiation field (Section 
3), quantification of the "order" of the atomic system (Section 4) and transfer 
of "order" from the radiation field to the atomic system (Section 5). 

The interesting point for solar surface magnetism is that weak magnetic 
fields (from 1 milligauss to 100 gauss, approximately) modify the atomic po- 
larization via the Hanle effect. However, in order to develop the Hanle effect 
as a reliable diagnostic tool of weak magnetic fields it is extremely important 
to fully understand the reference case of non-magnetic scattering polarization. 
With this motivation in mind Section 6 presents, for the first time, results of 
multilevel scattering polarization calculations taking fully into account all the 
relevant pumping mechanisms and the transfer of atomic polarization among all 
the levels involved. After considering in some detail the unmagnetized reference 
case a few examples in Section 7 will show how weak magnetic fields (from a few 
milligauss to a few gauss) modify the atomic polarization of multilevel atomic 
systems. Finally, Section 8 gives our concluding remarks. Section 2 is mainly 
dedicated to introducing the subject to those readers approaching the Hanle 
effect for the first time. 

2. Introduction to the Hanle effect 

This section is of introductory nature. Similar and additional information may 
be found in Hanle's (1924) paper on "Magnetic Effects on the Polarization of 
Resonance Fluorescence" , in the classical monograph of Mitchell and Zemansky 
(1934), in Landi Degl'Innocenti's (1992) contribution to the first IAC Winter 
School, in the book edited by Moruzzi and Strumia (1991), and in the pa- 
pers by Landolfi & Landi Degl'Innocenti (1986) and by Manso Sainz & Trujillo 
Bueno (1999). Other reviews and keynote articles related to this topic (includ- 
ing the determination of magnetic fields in solar prominences) can be found in 
Landi Degl'Innocenti (1990), Stenflo (1994), Faurobert-Scholl (1996) and Tru- 
jillo Bueno (1999). For a discussion concerning the possibilities of the Hanle 
effect for the detection of magnetic fields in stellar winds see Ignace et al. (1997). 

2.1. Hanle's "Doktorarbeit" 

The story began in 1923 when Wilhelm Hanle of Gbttingen University published 
the first correct interpretation of a previously observed phenomenon related to 
the effect of a weak magnetic field on the linear polarization of the spectral- 
line radiation scattered by mercury vapor illuminated anisotropically (see Hanle 
1923; 1924). With respect to the linear polarization corresponding to the zero 
magnetic field case, the observed influence of a weak magnetic field (of the order 
of 1 gauss) was a rotation of the plane of linear polarization (observed experi- 
mentally by Hanle himself) and a depolarization (clearly pointed out previously 
by Wood and Ellett in 1923). This so-called Hanle effect played a fundamental 
role in the development of quantum mechanics, since it led to the introduction 
and clarification of the concept of coherent superposition of pure states (see Bohr 
1924; Hanle 1924, 1925; Heisenberg 1925). The Hanle effect is directly related 



Atomic Polarization and the Hanle Effect 



3 



to the generation of coherent superposition of degenerate Zeeman sublevels of 
an atom (or molecule) by a light beamf]. As the Zeeman sublevels are split 
by the magnetic field, the degeneracy is lifted and the coherence (and, in gen- 
eral, also the population imbalances among the sublevels) are modified. As we 
shall illustrate below in the context of the solar polarized spectrum, this gives 
rise to a characteristic magnetic-field dependence of the linear polarization of 
the scattered light that is finding increasing application as a diagnostic tool for 
magnetic fields in astrophysics (see Astrophysical Spectropolarimetry , edited by 
Trujillo Bueno, Moreno Insertis and Sanchez 2001). 

2.2. The oscillator model for the Hanle effect 

Figure 1 is aimed at introducing the most basic ideas behind scattering polar- 
ization and the Hanle effect (see also Landi Degl'Innocenti, 1992). It is based 
on the classical interpretation of the Hanle effect for a triplet-type transition 
originating between a ground level with total angular momentum J/ = and an 
excited level with J u = 1. The atom is treated as a negative charge oscillating 
with angular frequency cuq and with a damping constant 7 given by the inverse 
of the lifetime of the excited atomic level (e.g. the upper level of a resonance 
line transition has a lifetime t\tf e ~l/A u [, A u \ being the Einstein coefficient for 
spontaneous emission). In solar-like atmospheres l/tiif e *C Ao;ii ne , with AuJn nc 
the width of the spectral line under consideration. 

As seen in Fig. 1, we are assuming that the atoms in the outer layers of 
the solar atmosphere are being illuminated by an unpolarized radiation field, 
which we are approximating by a unidirectional beam propagating in the radial 
direction. The unpolarized character of this radiation beam is indicated in Fig. 
1 by two perpendicular uncorrelated components of the electric field of the wave. 
Single-scattering events take place and the light polarization is measured for 
both forward scattering and 90 degree scattering as shown in Figure 1. 

Consider first the observation of the "north solar pole", where we have 
assumed in Fig. 1 that there is no magnetic field (i.e. there is no Lorentz force 
influencing the motion of the oscillating electron). Under these circumstances 
the atom can be represented by three independent linear oscillators vibrating at 
angular frequency loq along the axes of the reference system. As indicated in Fig. 
1, only the x and y oscillators are excited by the incident beam. The two excited 
oscillators radiate independently and decay radiatively with a damping constant 
7 = l/iiif e . If we observe along the direction of the incident beam (forward 
scattering case) we find that the radiation is obviously unpolarized. However, 
observing at 90° one finds that the radiation is linearly polarized along the x- 
axis, simply because the y-oscillator is seen pole-on. If we choose the positive 
reference direction for the Stokes Q parameter along the x-axis, we find Q = I 
and U = 0. Note that the same conclusion is obviously reached if the vibration 
of the x-oscillator is considered as that resulting from the coherent superposition 
of two counter-rotating circular oscillators, which are oscillating in phase with 
respect to each other at frequency cjq in the x—z plane. 



3 A coherent superposition of two or more sublevels of a degenerate atomic level is a quantum 
mechanical state given by a linear combination of pure states of the atomic Hamiltonian. 
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Figure 1. The oscillator model for the Hanle effect. 
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Consider now the observations of the "south solar pole", where we have 
assumed in Fig. 1 the presence of a weak magnetic field parallel to the solar 
surface and orientated along the y-axis. We have now to take into account the 
influence of the Lorentz force on the motion of the bound electron. The result 
is that the atom cannot be interpreted as three independent linear oscillators, 
but as a linear oscillator parallel to the magnetic field and two counter-rotating 
circular oscillators in the x—z plane oscillating at angular frequencies ujq + cjl 
and u>q — ujl, where wl = 8.79 x 10 6 i?(gauss) is the Larmor angular frequency 
(with B (gauss) indicating that the magnetic field is to be given in gauss). The 
resulting trajectory of the electron in the x—z plane is given by 



The trajectory described by these equations is an oscillation at frequency uq, 
with an amplitude that decays exponentially with a characteristic damping time 
given by tnf e = 1/7, and such that its oscillation axis is rotating at frequency 
• If the Zeeman splitting is sufficiently large so as to have wl > 1/ inf e (which 
can occur having still <C Au\i ne ) the oscillation axis can rotate several times 
before the oscillation amplitude is affected by the damping. The bound electron 
describes in the x—z plane the "daffodil" pattern shown in the lower part to 
the r.h.s. of Figure 1. Under these circumstances we will see totally unpolarized 
radiation for the 90 degree scattering case (i.e. for observation at the limb), but 
the maximum possible amount of linear polarization along the y-axis for forward 
scattering (i.e. for disc-centre observation). 

However, when the Zeeman splitting is similar to the natural width of the 
atomic level (i.e. when WL,~l/£iife) the oscillation axis rotates through an an- 
gle a within the characteristic damping time]]. The bound electron describes 
in the x—z plane the "rosette" shown in Figure 1. For the 90 degrees scatter- 
ing case (i.e. for observation at the limb) the y-oscillator is seen pole-on and, 
therefore, the observed polarization simply reflects the weighted average of the 
above-mentioned "rosette" pattern. With respect to the previous unmagnetized 
"north-pole case" , we now get (for observation at the limb) a depolarization and 
a rotation of the polarization plane (i.e. we now have a smaller Stokes Q value 
and a non-zero Stokes U signal). This rotation is counterclockwise for a magnetic 
field pointing toward the observer, as in Fig. 1, but clockwise if the magnetic 
field points in the opposite direction. Accordingly, for a mixed-polarity magnetic 
field topology within the spatio-temporal resolution element of the observations, 
the measured Stokes U parameter would be zero, while we would still be able to 
detect the depolarization effect by measuring Stokes Q. For a disc-centre obser- 
vation (forward scattering case) we would get (for the magnetic field orientation 
of Fig. 1) a net linear polarization signal along the y-axis, but of smaller am- 
plitude than that corresponding to the previous ^ 1 Anfe forward-scattering 
case. 



4 The rotation angle a = arctan([7/Q)/2 = ±arctan(2r)/2. For an atomic level of total angular 
momentum J the parameter T — gjuihtm e> with gj the Lande factor (which is unity for the 
upper level with J u = 1 of the assumed triplet-type transition). 





z 



(t) = ^e~ 7 * /2 sin(u; L t)cos(u;oi). 



(2) 
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Finally, it is also of interest to consider the "south pole" magnetized case, 
but assuming that the weak magnetic field is now orientated along the x-axis. 
It is easy to understand that all remains unchanged with regard to the forward 
scattering case (with the exception that now the observed linear polarization is 
along the x-axis, i.e. again along the magnetic field direction). However, for 
the limb observation we would only see the depolarization effect (i.e. a Stokes 
Q signal decreasing with the magnetic field, but U = always). Moreover, 
the Stokes Q signal does not vanish completely once the magnetic field becomes 
sufficiently large so as to have cjl 3> 1/tiife- Figure 2 corresponds to the situation 
of a magnetic field which is always orientated along the x-axis. Note that a 
magnetic field parallel to the solar surface and with an intensity in the saturated 
regime (i.e. wl 3> 1/iiife or T > 1) leads to a linear polarization amplitude 
for disc-centre observations that is an order of magnitude smaller than that 
corresponding to the unmagnetized reference case (r = 0) for an observation 
close to the solar limb (// = 0.1). Given the weakness of the "predicted" disc- 
centre polarization signals, two spectral lines of possible interest for a disc-centre 
Hanle-effect observational search of horizontal chromospheric fields would be the 
Ca I 4227 A line and the D 2 line of Ba II at 4554 Af]. 

2.3. The basic Hanle-efFect formula 

As we have seen, the Hanle effect produces a modification of the linear polariza- 
tion signals (quantified by the Stokes Q and U parameters) with respect to the 
reference case of non-magnetic scattering polarization. This occurs in a param- 
eter domain in which the transverse Zeeman effect is practically ineffective. The 
Hanle effect is sensitive to magnetic fields such that the corresponding Zeeman 
splitting is comparable to the inverse lifetime of the lower or the upper atomic 
levels involved in the line transition under consideration. The basic approxi- 
mate formula to estimate the maximum magnetic field intensity B (measured in 



where gj and t\a e are, respectively, the Lande factor and the lifetime of the 
atomic level under consideration (which can be either the upper or the lower 
level of the chosen spectral line transition). Depending on the astrophysical 
plasma under consideration, on the intensity and orientation of its magnetic 
field and on the spectral line chosen, we may find different Hanle-effect regimes. 
The most familiar one to astrophysicists in general is the "upper-level Hanle 
effect" , in which only the upper-level coherences are modified by the action of the 
magnetic field. However, we may also have a "lower-level Hanle effect" regime 
in which only the lower- level coherences are sensitive to the field. The most 



For having in principle the possibility of a positive detection the horizontal components of 
the chromospheric magnetic fields should not have a random azimuthal orientation within the 
spatio-temporal resolution element of the observation. 

6 The natural formula to write is uJLgj = 8.79xl0 6 B g j w l/tiif e . However, by omitting the 
8.79 factor we obtain an easier formula to remember (cf. Eq. 3), which informs us about the 
maximum magnetic field intensity to which the Hanle effect can be sensitive. 




10 6 B gj « 1/tufe , 



(3) 
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Figure 2. Hanle-effect radiative transfer simulation for limb and disc- 
centre observations. The figure shows the emergent fractional linear 
polarization versus the line frequency in units of the Doppler width. 
The positive reference-direction for the Stokes Q parameter is along 
the unit vector e2 of Figure 1. Assumptions: triplet- type transition 
of a two-level model atom with inelastic collisional destruction prob- 
ability e = 10 -4 in an isothermal model atmosphere and neglecting 
depolarizing elastic collisions. The magnetic field is parallel to the 
stellar surface and orientated along the x-axis of Figure 1. The sim- 
ulated observations are as in Fig. 1 (for fi = 1, but for /j, = 0.1 
instead of /x = 0). The intensity of the magnetic field is quantified by 
r = 8.79xl0 6 B {gauss) gj/A ul , with gj = 1. 
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general Hanle-effect takes place when the magnetic field affects the coherences 
of both levels simultaneously. Obviously, this dichotomy between an upper- 
level and a lower-level Hanle effect is only a suitable one if there exists a sizeable 
difference in the lifetimes of the lower and upper levels of the particular radiative 
transition. It certainly holds for solar spectral lines whose lower level is either 
the ground or a metastable level, as happens, for example, with the D-lines of 
Na I (Landi DegPInnocenti 1998) or with the Mg I b lines (Trujillo Bueno 1999). 
It is important to note that in the solar atmosphere (whose T c g~5800 K) the 
lifetime of a ground or metastable level is typically two orders of magnitude 
larger than the lifetime of the upper level. Therefore, for typical solar spectral 
lines the upper-level Hanle effect would be sensitive to magnetic fields of between 
1 and 100 gauss, while the lower-level Hanle effect could in principle be used for 
diagnosing much weaker fields, i.e. fields between 10 -3 and 1 gauss. Whether 
or not sub-gauss magnetic fields can actually exist in the highly conductive 
solar plasma is a question that most solar plasma physicists are inclined to 
answer negatively. The main result I aim at demonstrating in this article is that 
the "enigmatic" linear polarization signals observed by Stenflo et al. (2000) 
close to the solar limb in a number of chromospheric lines are undoubtedly 
due to the presence of atomic polarization in their metastable- lev el lower-levels. 
These observations of scattering polarization on the Sun have been confirmed 
recently (and extended to the Stokes U and V parameters) by Dittmann et al. 
(2001), Martinez Pillet et al. (2001) and by Trujillo Bueno et al. (2001) using 
different polarimeters attached to the Canary Islands solar telescopes. How 
such metastable-level atomic polarization can survive in the solar chromosphere 
is certainly a challenging question that I will also address in this paper. In 
particular, I will discuss briefly how a rigorous theoretical interpretation of this 
type of spectropolarimetric observations is giving us decisive new clues about the 
topology and intensity of the magnetic field of the "quiet" solar chromosphere. 

2.4. The utility of the Zeeman effect 

Before entering into details it is convenient to recall that the Zeeman effect is 
most sensitive in circular polarization (quantified by the Stokes V parameter), 
with a magnitude that scales with the ratio between the Zeeman splitting and the 
width of the spectral line, and in a way such that the V profile changes its sign 
for opposite orientations of the magnetic field vector. The longitudinal Zeeman 
effect in typical solar Fraunhofer lines is of course rather insensitive to sub-gauss 
magnetic fields. However, we should keep in mind that today's state-of-the-art 
polarimeters are perfectly able to detect Stokes V signals corresponding to a flux 
density of only a few gauss (see, e.g., Sanchez Almeida and Lites, 2000). Thus, 
unless this flux density is well below 1 gauss in the spatio-temporal resolution 
element of the observations, the longitudinal Zeeman effect itself should be con- 
sidered of complementary diagnostic value to the Hanle effect. In any case, it 
is important to emphasize that the Hanle effect, contrary to the Zeeman effect, 
does indeed work in any topologically complex weak-field scenario (i.e. even if 
the net magnetic flux turns out to be exactly zero) by producing a modification 
of the linear polarization signals (with respect to the zero magnetic field refer- 
ence case) that we can really "measure" if we succeed in rigorously modeling 
polarization signals due to multiple-scattering processes in the presence of weak 
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Figure 3. The Hanle effect versus the Zeeman effect. The figure 
shows the maximum polarization degree of a Ji = — > J u = 1 resonance 
line (with A u i = 10 7 s" 1 and Doppler width A^d = 10 10 s _1 ) versus 
the magnetic field intensity. The simulated observation is close to the 
limb (u, = 0.1) of the model atmosphere of Figure 2. The polarization 
signals are due either to the Zeeman effect (Z) or to the Hanle effect 
(H). The dashed-dotted line gives the sensitivity of the Hanle effect 
to a microturbulent and isotropic magnetic field. Note that for this 
kind of mixed-polarity scenario there is no Zeeman signal. The two 
other cases are for a magnetic field parallel to the stellar surface, and 
orientated either along the x-axis of Fig. 1 (curves with the symbol 
_L) or along the y-axis (curves with the symbol ||). The polarimetric 
signals corresponding to the transverse Zeeman effect (see the solid line 
labelled Z(_L)) have been multiplied by a factor 10 to make them visible 
in the figure. Note that for magnetic field intensities greater than 
10 gauss (i.e. for the saturated Hanle-effect regime where r^>l) the 
"Hanle-effect signal" is only sensitive to the orientation of the magnetic 
field vector, but not to its intensity. This can occur for a Zeeman 
splitting which is still a very small fraction of the width of the spectral 
line. (From Manso Sainz & Trujillo Bueno 2001b). 
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magnetic fields. In conclusion, besides being useful for "measuring" magnetic 
fields in solar prominences, the Hanle effect offers a promising diagnostic tool 
for the investigation of weak photospheric fields with mixed polarities over small 
spatial scales and for the exploration of chromospheric magnetic fields (because 
chromospheric lines are generally broad and the magnetic fields of the solar chro- 
mosphere relatively weak). Figure 3 illustrates what has just been pointed out 
in relation to the diagnostic interest of the Hanle and Zeeman effects. A suit- 
able illustration of the fact that the Hanle effect only operates in the Doppler 
core has been provided by Stenflo (1998) for the case of a Ji = 0^ J u = 1 line 
transition. 



3. How to quantify the order of the radiation field? 

In the standard multilevel radiative transfer problem (see Mihalas, 1978) the 
only quantity related to the radiative line transitions that plays a role in the 
statistical equilibrium equations is: 



7° - 

Jq — 



dx 



dtt 



(4) 



where x is the frequency measured from line centre in units of the Doppler 
width, 4> x the absorption line shape, O, the direction of propagation of the ray, 
and I the specific intensity of the radiation field, i.e. the Stokes / parameter. 
Regardless of the dependence of the radiation field on frequency and direction 
this quantity ( Jq) is always positive. Note also that Bi u Jq gives the contribution 
of the absorption process (between a lower level I and an upper level u) to the 
natural width of the lower level "Z" {B\ u is the Einstein coefficient for such an 
absorption process). If this lower level is either the ground or a metastable level, 
its lifetime is tnf e ~ 1/Bi u Jq if and only if the line transition I — > u is among the 
strongest ones. This applies, for instance, to the lower-levels of the Mg b lines, 
to the lower levels of the Ca n IR-triplet, and to many other spectral lines in the 
solar spectrum. All such lower levels are metastable. As a result, their lifetimes 
are about two orders of magnitude larger than the lifetime of the upper levels 
of the above-mentioned spectral lines. 

However, in the most general polarization transfer case, there are eight 
additional radiation field quantities that play a critical role in the statistical 
equilibrium equations (see Landi Degl'Innocenti, 1983). They are the spherical 
tensors (Jq ) of the radiation field (with K = 1, 2 and —K<Q<K) and they 
are the quantities we choose to quantify the "order" of the radiation field: 
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f r dCl [3 
4 = J dxf^—^-vV^, (8) 

where the azimuthal angle % an d H = cos# specify the orientation of each ray 
of direction fi. The reference direction for the Stokes Q and U parameters is 
situated in the plane perpendicular to and lies in the plane containing Q and 
the z-axis. The J-q components can be obtained easily from J^q = (— 1)®[Jq ]* 
(with Q > 0, and where the symbol "*" means complex conjugation.). 

In a weakly polarizing medium like the "quiet" solar atmosphere (cf. Sanchez 
Almeida and Trujillo Bueno 1999) these radiation field tensors are essentially 
related with the degree of "order" of Stokes /, with the degree of "order" of 
Stokes Q and U, and with the degree of "order" of Stokes V. By degree of 
"order" I mean, essentially degree of anisotropy, degree of breaking of the axial 
symmetry, and degree of "lack of antisymmetry" of the Stokes V parameter. An 
example of a completely "disordered" radiation field is that of a black body. It is 
isotropic, unpolarized and it has axial symmetry around any chosen direction in 
space. It is straightforward to show that its only non-zero radiation field tensor 
is the mean radiation intensity: Jq(v) = B u . However, the radiation field of 
a stellar atmosphere definitely has a degree of "order", which can be suitably 
quantified by the above-mentioned radiation field tensors. 

As an illustrative example, let us calculate the degree of anisotropy in an 
unmagnetized Milne-Eddington atmosphere (i.e. a medium in which the source 
function is S = a + br, with t the optical depth). Figure 4 shows, for various 
values of b/a, the variation with r of the anisotropy factor 

A = 4 (10) 

where Jq and Jq are given by Eqs. (4) and (5), respectively. To obtain the 
results of Fig. 4 one has to solve the transfer equation in order to calculate the 
variation with u. = cos9 of the Stokes I parameter {6 is the angle between the ray 
direction and the z-axis of the reference system, which we choose here along the 
normal to the stellar surface) . Note that A < for the case of an atmosphere with 
no gradient in the source function (i.e. the curve with b/a = corresponding to 
S = a). The anisotropy factor {A) is essentially negative in atmospheres with 
very small b/a values (see, for example, the curve with b/a = 0.1), while *4>0 
in atmospheres with sufficiently large b/a ratios (see, for example, the curves 
with b/a>l). Note also that the larger b/a the larger the anisotropy factor. 

These results can be understood intuitively as follows. First, note that the 
tensor Jq (cf. Eq. 5) is dominated by the first term having the Stokes / contri- 
bution. Therefore, predominantly vertical rays (i.e. those with fi > l/\/3) make 
positive contributions to Jq, while predominantly horizontal rays (i.e. those with 
/i < l/\/3) make negative contributions^. Second, as illustrated in Fig. 5, in a 



7 The angle corresponding to fi = is known as Van Vleck's angle. 
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Figure 4. The anisotropy factor in Milne-Eddington atmospheres. 
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Figure 5. The anisotropic illumination in a stellar atmosphere. 
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stellar atmosphere the outgoing intensities diminish with decreasing n. (i.e. they 
show limb darkening), while the incoming intensities augment with decreasing 
/j, (i.e. they show limb brightening). In other words, in a stellar atmosphere 
the outgoing radiation is predominantly vertical, while the incoming radiation is 
predominantly horizontal. Thus, the outgoing intensities tend to produce posi- 
tive contributions to Jq , while the incoming intensities tend to produce negative 
contributions. Therefore, there is competition. It wins the subset of intensities 
(outgoing or incoming) having the largest variation with \x. If b/a = the source 
function is constant (S = a) and the outgoing intensities have no \x dependence. 
This implies that .A<0 because the incoming intensities always show a depen- 
dence with /i, even for a constant-property atmosphere. This is simply due to 
the presence of the stellar surface. However, as soon as b/a becomes "large 
enough" the limb darkening of the outgoing intensities becomes more important 
than the limb brightening of the incoming intensities and the anisotropy factor 
(A) becomes positive. It is straightforward to show that, for a given temperature 
gradient of the model atmosphere, the corresponding source function gradient 
is larger the bluer we go in the spectrum. This is the main reason why the 
anisotropy factor generally increases as we move in the spectrum from the IR 
toward the UV. 

The anisotropy factor (cf. Eq. 10) is a fundamental quantity in scattering 
polarization. It is also of interest to note that its possible values are bounded as 
dictated by the following expression: 

i 1 

<A<^. (11) 



2^/2 ~ ~ y/2' 

In the chosen reference frame (with the z-axis along the radial stellar direction), 
the largest A value corresponds to an illumination coming from purely vertical 
(radial) radiation beams and the lowest one to a purely horizontal radiation 
field. 

The remaining radiation field tensors are more subtle, but their physical 
meaning can also be understood intuitively. The tensors Jq (with Q = 1 and 
Q = 2) are complex quantities. They quantify the breaking of the axial sym- 
metry of the Stokes /, Q and U parameters. If the physical properties of the 
stellar atmosphere model depend only on the radial direction, the ensuing Jq 
tensors (with Q = 1 and Q = 2) are zero unless the magnetic field is inclined 
with respect to the radial direction. Figure 6 shows the variation with depth 
in a stellar model atmosphere of the real and imaginary parts of the Jq tensors 

normalized to Jq. Note that Jq / Jq <C 1, i.e. that the degree of "order" of the 
radiation field which interacts with the atoms of the stellar atmospheric model 
is weak. This must imply that the degree of "order" of the atomic system (i.e. 
the atomic polarization) has also to be weak. These considerations can suitably 
be taken into account in order to facilitate some analytical insight into and the 
numerical solution of the general equations. 

The tensors Jq are related to the circular polarization. Note that Jq = 
if the Stokes V profile is perfectly antisymmetric, and that J\ = if there is no 
azimuthal % dependence in the Stokes V parameter. We thus need to have net 
circular polarization in order to permit non-zero values for the Jq tensors. 
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Figure 6. The tensors Jq/Jq corresponding to the self-consistent so- 
lution in the same stellar atmospheric model of Figure 2. The assumed 
horizontal magnetic field has T = 1 and is orientated at 30° with re- 
spect to the x-axis of Figure 1. The solid lines give the exact result. 
The dashed line of the panel for Jq / Jq neglects the Stokes Q contribu- 
tion of Eq. (5). The dashed lines of the remaining panels neglect the 
contribution of the Stokes I parameter to Eqs. (6) and (7). Therefore, 
the Stokes I parameter is essential for calculating Jq, while the Stokes 
Q and U parameters cannot be neglected (in general) when calculating 
J\ and Jf. 
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The most general situation in which all the above-mentioned radiation field 
tensors are non-zero is the case of a magnetized stellar atmosphere with macro- 
scopic velocity gradients. In this article I will consider static stellar atmospheres 
for which the Jq tensors can be assumed to be zero. As shown in Section 5, I 
will focus mainly on the generation of atomic alignment, which is directly related 
to the Jq tensors and with the Hanle effect (i.e. with the modification of the 
linear polarization signals due to weak magnetic fields) . In a future occasion I 
will address the issue of the generation of atomic orientation, which is intimately 
related to the Jq tensors and with the modification of the circular polarization 
signals. 

4. How to quantify the "order" of the atomic system? 

The "order" of the atomic system (i.e. the atomic polarization) can be suitably 
quantified via the atomic density operator (p A ) of quantum mechanics. This 
operator is represented in the basis {|n}} of energy eigenvectors via a matrix 
called the atomic density matrix whose elements are: 

As clarified below, the diagonal elements (p^) quantify the population of the 
state \i), while the non-diagonal elements (p{j, with i ^ j) quantify the degree 
of quantum interference (or coherence) between the states \i) and \j). 

Let us give some information regarding the exact meaning of the atomic 
density matrix. First, note that the physical system under consideration is 
a volume element of a stellar atmosphere, which can be considered as being 
composed of three subsystems: the atoms (A) which emit, absorb, or scatter the 
radiation, the material perturbers (P) capable of influencing the excitation state 
of the atoms through inelastic and elastic collisions, and the radiation field (R) 
itself. The quantum-mechanical Hamiltonian operator H of the whole system 
can be written, in the Schrbdinger representation, in the form H = Ho + V, 
where Ho = Ha + Hr+Hp is the unperturbed Hamiltonian made up of the sum 
of the energies of the atoms, radiation and perturbers, while V accounts for the 
various possible interactions among atoms, perturbers and radiation. Second, 
recall that there are two reasons which lead to the introduction of probabilities 
in the description of the state of a physical system such as that just defined. 
One is the quantum-mechanical uncertainty related to the measurement process. 
The other is due (as in classical statistical mechanics) to our lack of complete 
information on the initial state of the system. This means that its description 
cannot be given in terms of a pure state |*), but through a statistical mixture of 
states (with the probability p^ of finding the system in the state vector \^k))- To 
incorporate into the quantum mechanical formalism the incomplete information 
we posses about the state of the system, the mixed state density operator (i.e. 
p = J2kPk\^k)(^k\) was introduced (von Neumann, 1927; but see Fano, 1957). 
The essential point is that the microscopic interactions of our atomic subsystem 
with the perturbers and the radiation field drive the atomic subsystem into a 
mixed state which can be described by the density operator but not by a pure 
state. 
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To grasp the usefulness of the atomic density operator, consider a complete 
set {\n)} of eigenvectors of the atomic Hamiltonian (i.e. H^n) = E n \n), with 
E n the energy eigenvalue corresponding to |n)). To be able to make physical 
predictions about measurements bearing only on the atomic subsystem (A), the 
density operator p A of this subsystem is introduced (i.e. p A = Trp (Trpp), 
where the symbols Trp and Trp mean the traces over the perturbers (P) and 
the radiation field (R) coordinates respectively; note that the trace of a matrix 
is the sum of its diagonal elements). Calculating the matrix elements p A - of p A 
in the set {|n)} one finds: 

• p A represents the average probability of finding the system in the state | i) 
(which implies that p A quantifies the population of the state 

• p A j (with i ^ j) gives account of the interference effects between the states 
\i) and \j), which can appear simply because each state > of the 
statistical mixture is (in general) given by a coherent linear superposition 
of the states As mentioned above, the non-diagonal elements of p A 
are called coherences. 



The density operator p of the total system (atoms, perturbers and radia- 
tion) contains all physically significant information on the system. In fact, the 
expectation value of any observable, described in the Schrodinger picture by a 
Hermitian operator O, is given by (cf. Fano, 1957): 



(0)(t) = Tr(pO) = Tr(p J / ), (13) 

where O/ and pi are the corresponding operators in the interaction (or Heisen- 
berg) picture of quantum mechanics (see any suitable advanced textbook). As 
indicated by this expression, the expectation value of an observable has the 
same structural form in the two quantum mechanical approaches (i.e. in the 
Schrodinger and in the Heisenberg representations). However, the Heisenberg 
picture has the advantage of forcing the time dependence of wave functions to 
arise solely from the effect of the perturbing Hamiltonian V, thus facilitating a 
perturbation approach. 

Equation (13) is one of the two basic equations on which the density-matrix 
polarization transfer theory is based. The other fundamental equation is the 
equation of motion for the density operator pj, which reads^j 

The combination of Eqs. (13) and (14) leads to an exact equation for 
d{0)/dt, which can be used to derive the time evolution of physical quantities 
(e.g. the time evolution of the diagonal and non-diagonal elements of the atomic 
density matrix p A due to the coupling of the atomic system A with the radiation 



It is often called the Liouville equation in the interaction picture. The corresponding equation 
in the Schrodinger picture is identical to Eq. (14), but with p instead of pi and with the total 
Hamiltonian H in place of Vi. It can be easily derived from Schrodinger's equation. 
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field R considered as a reservoir; cf. Cohen- Tannoudji, 1977). As reviewed in 
more detail by Trujillo Bueno (1990), the above-mentioned exact equation for the 
time evolution of the observables is the starting point in Landi Degl'Innocenti's 
(1983) QED derivation of the statistical equilibrium and Stokes-vector transfer 
equations. 

Finally, we should mention that there are two relevant representations of 
the atomic density operator: the standard (already introduced above) and the 
spherical statistical tensor representation (Omont, 1977; Blum, 1981). Let us 
introduce them in terms of the basis of eigenvectors of the angular momentum 
(\aJM >, with a indicating the quantum numbers of a given term — e.g. a = 
nLS if the atom is described by the L — S coupling scheme). 

In the standard representation the coherences between magnetic sublevels 
pertaining to the same J-level are given by[]: 

p%j(M,M') = (aJM\p A \aJM'}, (15) 
while the total population of the J-level is 



nj 



J2PaAM,M), (16) 



M 



where p^j(M, M) is the population of the sublevel with magnetic quantum num- 
ber M. 

In the spherical statistical tensor representation the density-matrix elements 
are denoted by the symbol pq (with K = 0, ...,2J and —K<Q<K). The Pq 
elements are given by the following linear combinations of the density-matrix 
elements of the standard representation (cf. Omont, 1977): 

pq^J) = E (-i) J - M V2KTT ( 3 _ J , K ) pUmM), (17) 

MM' V W 1 

where the 3-j symbol is defined as indicated by any suitable textbook on quantum 
mechanics. 

For instance, for a level with total angular momentum J = 1 

Po = ^IPii 1 ' !) + ^(°. °) + Pii' 1 ' = 7^ n J=^ ( 18 ) 
pj = -^(l,l)-pf(-l,-l)], (19) 

pl = --^[^(1,0) +^(0,-1)], (20) 

pI = ^br(i, i) - 2pr(o, o) + ^(-1, -i)], (21) 



9 We may also have coherences between sublevels of different J levels. In some cases (e.g. the H 
and K lines of Ca n) such coherences may have some observable effects (see Stenflo, 1980). 
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combinations of the populations of the various Zeeman sublevels corresponding 
to the level of total angular momentum J. The total population of the atomic 
level is quantified by \JU + Tp(], while the population imbalances among such 
Zeeman sublevels are quantified by p\ (i.e. by the Q = orientation coefficient) 
and by p\ (i.e. by the Q = alignment coefficient). However, the Pq elements 
with Qt^O are complex numbers given by linear combinations of the coherences 
between Zeeman sublevels whose magnetic quantum numbers differ by Q. In 
fact, since the density operator is Hermitian, we have that for each spherical 
statistical tensor component pg with Q > 0, there exists another component 

with Q < given by p^g = (— 1)^[/9q]*- We thus have (2J+ l) 2 density-matrix 
elements corresponding to each level of total angular momentum J (both in the 
standard and spherical tensor representations). 

Both formalisms are totally equivalent. However, in the context of scatter- 
ing polarization and the Hanle effect, it is more convenient to work within the 
framework of the spherical statistical tensor representation due to the following 
reasons (cf. Landi Degl'Innocenti 1982). The first advantage is directly related 
to the fact that the density matrix elements on the basis of the eigenvectors of 
the angular momentum depend on the reference system chosen to define such 
eigenvectors. It turns out that the spherical statistical tensors defined by Eq. 
(17) have simpler transformation laws with respect to rotations of the coordinate 
system: their transformation law involves just one rotation matrix instead of the 
product of two rotation matrices. The second advantage is that the pg elements 
themselves provide the most suitable way of quantifying, at the atomic level, the 
information that we need to be able to calculate all the "sources" and "sinks" of 
linear and circular polarization within the medium under consideration (i.e. they 
have a clear physical interpretation). For instance, if /?q( J M )^0 and Pq(J u )^0 we 
have local sources of linear and circular polarization, respectively, even in the 
absence of magnetic fields. In the Hanle effect regimes considered in this article 
the width of the spectral line is much larger than the Zeeman splitting and the 
Stokes Q and U components of the line emission vector are given by: 



eg = €Q W { j ) uJi ^^-j={p 2 - l)pl - V3p^Jl - p 2 (cosxRe[p?] - sin^Im^ 2 ]) 

- ^-(1 + p 2 )(cos2 X Re[p 2 ] - sin2 X Im[/> 2 ])},(24) 



e v = e w£ } J; VS { \jl - p 2 (sin X Rc {p\\ + cosxIm[p 2 ]) 




(25) 
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where the Pq values are those of the upper level of the line transition under con- 
sideration, eo = (hv/4iT)A u i4> x J\f\/2J u + 1 (with J\f the total number of atoms 

(2) 

per unit volume), Wj^j is the symbol introduced by Landi Degl'Innocenti (1984) 
(which depends only on J u and J;), and where the orientation of the ray is spec- 
ified by \i = cos6* (with 9 the polar angle) and by the azimuthal angle %. The 
elements r)Q and r\\j of the absorption matrix are given by identical expressions 
(i.e. by r] Q = e Q and by rftj = e v ), but with r/o = {hv j Am)B lu $ x M + 1 

instead of eo, wffj u instead of w^ Ji and with the pq values of the lower level 
of the line transition (instead of those of the upper level). Note that cq and tjq 
depend on both the population imbalances (pg) an d on the coherences (pQ, with 
Q = 1, 2), while ejj and r/;y depend only on the coherences. Recall also that the 
emission vector and the absorption matrix are the two fundamental quantities 
which appear in the Stokes-vector transfer equation. 

Finally, an additional reason for the suitability of the spherical tensor rep- 
resentation is that the limiting case in which polarization phenomena are ne- 
glected (cf. Mihalas, 1978) can be obtained simply by retaining only the terms 
with K = Q = 0. 

5. Transfer of "order" from the radiation field to the atomic system 

There are two mechanisms which can lead to transfer of "order" from the ra- 
diation field to the atomic system: upper-level selective population pumping 
and lower-level selective depopulation pumping. The requirements are that the 
pumping light must be necessarily anisotropic, and/or polarized and/or to have 
spectral structure over a frequency interval Av smaller than the frequency sep- 
aration between the Zeeman sublevels. 

Upper-level population pumping occurs when some upper-state sub- 
levels have more chances of being populated than others. For instance, if an 
unpolarized light beam propagating along the direction chosen as the quantiza- 
tion axis illuminates a gas of two level atoms with Ji = and J u = 1, only the 
transitions corresponding to AM = ±1 are effective, so that no transitions occur 
to the M = sublevel of the upper level. Thus, in the absence of any relaxation 
mechanisms, the upper- level sublevels with M = 1 and M = — 1 would be more 
populated than the M = sublevel and the fractional upper-level alignment 
coefficient (Tq(u) = Pq(u)/pq(u) = l/y/2. 

Lower-level depopulation pumping occurs when some lower-state sub- 
levels absorb light more strongly than others. As a result an excess population 
tends to build up in the weakly absorbing sublevels. For instance, if an unpo- 
larized light beam propagating along the direction chosen as the quantization 
axis illuminates a gas of two level atoms with J; = 1 and J u = 0, only the 
transitions corresponding to AM = ±1 are effective, so that no transitions can 
occur out of the M = sublevel of the lower level. On the other hand, the 
spontaneous de-excitation from the upper level populates with equal probability 
the three sublevels (M = — 1,0, +1) of the lower level. In the absence of any 
relaxation mechanisms, the final result of this optical-pumping cycle is that all 
atoms will eventually be pumped into the M = sublevel of the lower level, 
and the medium will become transparent (Happer, 1972). Under such ideal 
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laboratory conditions (illumination by a unidirectional light beam and absence 
of depolarizing mechanisms) the fractional lower-level alignment coefficient is 

°o(0 = Po(0/Po(0 = ~~ v^, i-e. a factor 2 larger (in absolute value) than the 
fractional upper-level alignment corresponding to the previous triplet-line case. 

The Nobel laureate Alfred Kastler (1950) was actually the first scientist 
to propose that optical pumping under laboratory conditions can be used as a 
method to change the relative populations of Zeeman sublevels and of hyperfine 
levels of the ground state of atoms. What we have been emphasizing over the last 
few years is that the very same mechanism (lower-level depopulation pumping) 
is operating in the atmospheres of the stars, and that it constitutes an essential 
physical ingredient for understanding the second solar spectrum. 

Lower-level depopulation pumping in solar-like atmospheres was investi- 
gated in detail by Trujillo Bueno and Landi DegFInnocenti (1997) including 
depolarizing collisions and radiative transfer effects. We chose that particular 
type of line transition (i.e. 1 — > — > 1) to point out clearly that lower-level 
depopulation pumping due to the anisotropy of the solar radiation field can lead 
to sizeable amounts of ground-level atomic polarization and to emergent linear 
polarization signals with amplitudes in the observable range. Soon afterwards, 
Landi Degl'Innocenti (1998) succeeded in deriving analytical expressions for the 
Stokes-Q component of the emission vector and of the absorption matrix cor- 
responding to the Di and D2 lines of Na I, and could show by adjusting free 
parameters that a certain amount of atomic polarization in the hyperfine com- 
ponents of the ground level of sodium leads to a remarkably good fit of the 
complex fractional linear polarization pattern observed by Stenflo and Keller 
(1997), including the "enigmatic" line-centre peaks. A subsequent theoretical 
investigation (Trujillo Bueno 1999) was aimed at demonstrating that depop- 
ulation pumping in the solar atmosphere (and the induced lower-level atomic 
polarization) can actually produce similar amounts of linear polarization for 
some groups of line transitions (having different J\ and J u values) for which 
the simplified resonance line polarization theory (which neglects the influence of 
lower-level depopulation pumping) predicts drastically different emergent polar- 
izations. Based on this result, I suggested a possible explanation of the "enig- 
matic" linear polarization amplitudes observed by Stenflo et al. (1983, 2000) in 
the Mg I b lines (see Trujillo Bueno 1999; Sections 6 and 7). That theoretical 
investigation demonstrated that the presence of a sizeable amount of atomic po- 
larization in the metastable lower-levels of the Mg b\ and 62 lines would explain 
in a natural way the similar polarization amplitudes observed in the Mg b lines 
(and also those observed in the three lines of multiplet No. 3 of Ca I, which 
is also a 3 P° — 3 S multiplet). This result for the Mg b lines was particularly 
encouraging because, contrary to the case of the sodium D-lines, it is totally 
impossible to argue that there might perhaps exist an alternative explanation of 
the observed polarization peaks based on a multilevel scenario characterized by 
the absence of lower-level polarization (see Trujillo Bueno 1999; Sections 5 and 
6). 

In addition to the two previous pumping mechanisms, which allow the di- 
rect transfer of "order" from the radiation field to the atomic system, there is an 
additional pumping process which also comes into play: repopulation pumping. 
This pumping occurs either when the lower-level is repopulated as a result of 
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spontaneous decay of a polarized upper level or when the upper-level is repopu- 
lated as a result of absorptions of a polarized lower level. As can be intuitively 
expected, the repopulation pumping rates are proportional to the atomic polar- 
ization of the atomic levels. 

5.1. The two- level atom in the absence of magnetic fields 

The best thing we can do to facilitate understanding of all these pumping mech- 
anisms is to show them in action for a particularly illuminating example: the 
scattering line polarization problem in a one-dimensional unmagnetized stellar 
atmosphere formed by a gas of two level atoms with J\ = J u = 1 (cf. Trujillo 
Bueno 1999). For this case, the only non-zero radiation field tensors are Jq 
and Jq and the unknowns of the problem at each spatial point are simply Pq(1), 
Po(0> Po( u ) an d Po( u )- The rate equations which govern the temporal evolution 
of the alignment coefficients (pg) of the upper and lower levels can be deduced 
applying the density matrix polarization transfer theory of Landi Dcgl'Innocenti 
(1983). We obtain: 



Tt pl{u) 



B, 



B, 



lu jo o 



B. 



lu 72 2 / 



V2 



-A ulP l(u) (26) 



+C^p 2 (l)-(C ul + D u )pl(u), 



dt 



pl(l) = -A uW l{u) + 



Bi u JqPq{1) - 2Bi u J p (l) 



B, 



lu 72 2 



y/2 



(27) 



+2C ulP l{u)-2(C lu + D l )pl(l). 



These equations have contributions coming from transfer and relaxation 
rates due to radiative and collisional terms. I have neglected the stimulated 
emission terms. Assuming statistical equilibrium (i.e. (ipq/dt = 0), neglecting 
also inelastic collisions (i.e. the C/ u and C u \ terms) and the terms containing the 
products JqPq (because in a weakly polarizing medium like the solar atmosphere 
J 2 /J ° < 1 and pi/pi < 1) we find: 
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2 (1 + 5u) 
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Po(low) 



Iu-Jq 



(28) 



(29) 



where 5 U = D u /A u i and 5i = Di/Bi u Jq, with D u and Di the depolarizing rates 
due to elastic collisions. 

The first term within the brackets of Eq. (28) is due to the upper-level 
population pumping mechanism, while the first term within the brackets of Eq. 
(29) is due to the lower-level depopulation pumping mechanism. Both are given 
by the rate B\ u JqPq(Iow); the only difference is that the upper-level pumping 
contribution is multiplied by the lifetime of the upper level (i.e. by 1/A u i), while 



22 



Javier Trujillo Bueno 



the lower-level pumping contribution is multiplied by the lifetime of the ground 
level (i.e. by 1/Bi u Jq, which for optical line transitions in solar-like atmospheres 
is about two orders of magnitude larger than the upper-level lifetime) . 

The second term within the brackets of Eq. (29) is due to the repopulation 
pumping resulting from the spontaneous decay of the polarized upper level, while 
the second term within the brackets of Eq. (28) is due to the repopulation 
pumping resulting from the absorption process from the polarized lower level. 

The previous two equations have been useful in identifying the contribu- 
tions corresponding to the various pumping mechanisms. However, the relevant 
quantities of interest to understanding the observed polarization amplitudes are 
the fractional atomic polarization (<7q = Pq/Po) °f each level. In fact, the emer- 
gent fractional linear polarization close to the solar limb is approximately given 
by: 

j « Wag (up) - .Sag (law), (30) 

where W and Z are simply numbers given by the total angular momentum values 
of the lower and upper atomic levels of the line transition under consideration^]. 
In this formula the Ug values have to be calculated at the spatial point (situated 
along the line of sight) where the optical depth is unity. This formula (see 
Trujillo Bueno 1999) may be considered as the generalization of the Eddington- 
Barbier relation to the non-magnetic scattering polarization case. (Note that 
I have rewritten it here in a way such that the positive reference-direction for 
the Stokes Q parameter is now chosen along the line perpendicular to the radial 
direction through the observed point, and not along the radial direction as it 
was chosen in the above-mentioned paper and in Eqs. 24 and 25.) 

The equations for Oq(up) and 0q(1ow) can be obtained easily after divid- 
ing Eqs. (28) and (29) by p®(l) and taking into account that Pq(u)/po(1) ~ 
Bi u J$/Aui. (Note that this is a good approximation for a weakly anisotropic 
medium in which the inelastic collisional rates have been assumed to be much 
smaller than the radiative rates.) The result reads: 



ag(up): 



0q (low) 



-1 



2 (1 + 5 U ) 



2(1 + S t ) 



-^ + 0-0 (low) 



|--o 2 (up) 



(31) 



(32) 



We thus see that there exists a fascinating closed loop connecting the lower- 
level and upper-level polarizations. The atomic polarization of the upper level 
produced by the anisotropic illumination of the atoms is modified (typically 
enhanced!) because of the atomic polarization of the lower level, and vice versa. 
The only way to frustrate this remarkable communication between the two levels 
is by forcing one of the two levels to be totally unpolarized, which can only occur 
in practice if either 5 U or 5i turns out to be very much larger than unity. In fact, 



10 Note that W = Wj^j, and Z = Wj^j , where the symbol Wj, is the one introduced by Landi 
Degl'Innocenti (1984). For instance, for a 1— +1 transition W = Z = —1/2. 
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if we assume that the atomic polarization of the lower level is totally destroyed 
by elastic collisions (i.e. 5[ — > oo) and that the upper level is insensitive to such 
collisions (i.e. 5 U = 0) we find <7q(up) rj — 0.5(Jq/Jq), in agreement with the 
self-consistent numerical results given by the dashed lines of Fig. 1 of Trujillo 
Bueno (1999). However, if we assume that both atomic levels are insensitive to 
elastic collisions (i.e. that Si = S u = 0), we then find ctq (up) — <7q(1ow) rj 
— Jq/Jq, in agreement with the self-consistent numerical result given in Fig. 3 
of Trujillo Bueno (1999). 

5.2. The two-level atom in the presence of weak magnetic fields 

What do we obtain for the atomic polarization of the ground and excited levels 
of a 1— >1 transition if we include the effect of a weak magnetic field? Obviously, 
the situation becomes more complicated. First, the number of unknowns at each 
spatial point is 12 instead of just four. The density matrix elements whose values 
are sought aref[3 Po(l), Po(l), the real and imaginary parts of pf(l) and of 
Pq(u), Pq(u) and the real and imaginary parts of p\(u) and of p^u). Second, 
there are two "natural" choices for the direction of the z-axis of the reference 
system used to formulate the rate equations: (a) the magnetic field reference 
frame in which the z-axis is aligned with the local magnetic field direction and (b) 
the stellar radial direction itself. For developing a general multilevel scattering 
polarization and Hanle-effect code we prefer option b (which is in fact the choice 
in Sections 6 and 7). However, as done in the remaining part of this section, 
it is convenient to formulate the rate equations in the magnetic field reference 
frame if the interest lies in achieving some analytical insight. 

The two following expressions are a particular case of the equations derived 
by Landi DegFInnocenti (1985). They are valid for a 1— >1 transition of a two- 
level atom in a weakly anisotropic and magnetized stellar atmosphere in which 
the effect of collisions is assumed to be negligible: 

■*<"> " " - ( i + K' + °£gj-i M)C f ' (33) 

These formulae are valid for K = 1 and K = 2. They show that, once the 
fractional radiation field tensors (i.e. J_n/ J$) are known in the magnetic field 
reference frame, the fractional atomic polarization of each level depend only on 
the dimensionless quantities T u and Ti, which are given by: 

T u = 8.79 x W e BgjJA ul , (35) 
= 8.79 x ltfVgjjB lu Jl (36) 



1 Actually, there are six additional unknowns which are related with the atomic orientation (i.e. 
with the Pq's), but these density-matrix elements are zero if the radiation field tensors Jq = 0. 
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The previous equations for o~q (u) and <Tq (I) show that in the magnetic field 

reference frame the population imbalances (i.e. a^) are insensitive to the mag- 
netic field, while the coherences (i.e. <7q with Qj^O) are reduced and dephased 
as the magnetic field is increased. Finally, note that magnetic fields such that 
~ 1 are very efficient in depolarizing the ground level, which has a significant 
feedback on the upper-level coherences, while fields such that T u ~ 1 are very 
efficient in depolarizing the excited level. 

Some polarization diagrams illustrating the lower-level and upper-level Hanle 
effects for single scattering events in two-level atomic models have been pre- 
sented by Landolfi & Landi Degl'Innocenti (1986). Similar type of Hanle-effect 
diagrams, but taking into account radiative transfer effects in isothermal stellar 
atmospheres, will be published by Manso Sainz & Trujillo Bueno (2001b). 

6. Multilevel scattering polarization 

The theoretical interpretation of scattering polarization signals requires a calcu- 
lation of the the atomic polarization considering complex atomic and molecular 
systems. This is a very involved non-linear and non-local problem which has 
been solved with sufficient generality only recently (see below). It consists in 
calculating, at each spatial grid-point of the assumed astrophysical plasma, the 
density-matrix elements which are consistent with the properties of the polarized 
radiation field generated within the medium. 

From multilevel radiative transfer simulations without polarization physics 
we know that the two-level atom model is, with few exceptions, an unsuit- 
able approximation for modeling the Fraunhofer spectrum. Obviously, the same 
statement should be applied (but with much more emphasis) to the realm of the 
second solar spectrum, which is due to coherence transfer effects. In order to 
exploit its rich diagnostic potential (e.g. with the aim of exploring magnetotur- 
bulence and/or to improve our understanding of chromospheric magnetism) it 
is crucial to investigate the generation and transfer of atomic polarization con- 
sidering realistic multilevel models and taking fully into account all the above- 
mentioned pumping mechanisms in the presence (and absence) of the Hanle 
effect induced by weak magnetic fields (e.g. from 1 milligauss to 100 gauss). 
During the last year Rafael Manso Sainzj^] and I have successfully developed a 
general multilevel scattering polarization radiative transfer code, which is cur- 
rently allowing us to carry out this type of investigations via realistic numerical 
simulations. Our multilevel scattering polarization code solves the following set 
of equations: 

(a) The rate equations giving the time evolution of the atomic density 
matrix elements. In these equations the radiative rates are given in terms of 
the radiation field tensors introduced in Section 3. For the moment we assume 
statistical equilibrium (SE) (i.e. dpq/dt = 0). 

(b) The Stokes- vector transfer equations where the components of the emis- 
sion vector and those of the absorption matrix depend on the density-matrix 
elements themselves. For the moment, we assume one-dimensional (ID) plane- 



He is now at the University of Florence, after having completed his Ph.D. thesis at the IAC. 
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parallel geometry, but we have also developed suitable formal solution methods 
for considering spherical circumstellar envelopes and realistic 3D model atmo- 
spheres with horizontal inhomogeneities. 

These SE and RT equations were derived from the principles of quantum 
electrodynamics by Landi Degl'Innocenti (1983). This density-matrix polariza- 
tion transfer theory is based on the Markovian assumption of complete frequency 
redistribution (see a critical review of this CRD polarization transfer theory in 
Trujillo Bueno 1990). For problems where the only significant coherences are 
those between the sublevels of degenerate levels, this CRD theory should provide 
a physically consistent description of scattering phenomena if the spectrum of 
the pumping radiation is flat across a frequency range wider than the inverse 
lifetime of the levels (Landi Degl'Innocenti et al. 1997). As we shall see below, 
this seems to be a sufficiently good approximation for modeling the observed 
polarization in a number of spectral lines of diagnostic interest (e.g. the Ca n 
IR triplet and the Mg b lines)0. 

Other scientists (e.g. Bommier and Sahal-Brechot 1978) had previously for- 
mulated the statistical equilibrium equations for the Hanle effect regime applying 
also the theory of the master equation of a "small system" (the atom) interacting 
with a "big reservoir" (the radiation field) (cf. Cohen- Tannoudji 1977; Cohen- 
Tannoudji et al. 1992). The interest of Landi's theoretical work stems from 
the fact that he was able to achieve a general unified theoretical derivation of 
both, the statistical equilibrium equations and the radiative transfer equations^. 
From the point of view of numerical radiative transfer our coherency transfer 
code is based on the fast iterative methods which the author of this article de- 
veloped for multilevel scattering polarization and Hanle effect applications (see 
Trujillo Bueno 1999; and note that the two-level atom scattering polarization 
problem with lower-level depopulation pumping solved in that paper is also a 
non-linear and non-local problem of the same nature as the general multilevel 
case considered here). 

6.1. Three-level systems and the "enigmatic" Ca n 8662 A line 

The quest to understand the physical origin of the second solar spectrum will 
be facilitated if we first gain some physical insight on how the various pumping 
mechanisms leads to atomic polarization in three-level atomic systems under 
stellar atmospheric conditions. 

There are a number of schematic three-state configurations which our laser 
physics colleagues like to include in their textbooks: cascade, vee (V) and 
lambda (A). The linkage pattern is 1 <-» 2 <-> 3 for any ranking of the rel- 
ative energy values. Thus, the cascade configuration is the one in which the 
energies increase with index i, E\ < E% < E%. The vee configuration occurs 



In any case, it cannot be emphasized sufficiently the importance of continuing the investigations 
to generalize the density-matrix polarization-transfer theory to partial frequency redistribution 
and to situations where the excitation of the atomic system is sensitive to the spectral structure 
of the pumping radiation (see Bommier 1997 a,b; Landi Degl'Innocenti et al. 1997). 

It is of historical interest to mention the works of House (1970 a,b; 1971), who was among the 
first to provide some significant contributions to the theory of scattering polarization and the 
Hanle effect. See also the work of Litvak (1975) on polarization of astronomical masers. 
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Figure 7. Left: The fractional atomic polarization of the lower level 
of the 8662 A line versus its optical depth (see the thin solid line 
yielding a surface value of —3%). The thick solid line yielding a surface 
value of —7% shows — Jq/Jq. Right: The corresponding emergent 
fractional linear polarization of the 8662 A line at ^ = 0.1. 



when the central state (i.e. level 2) lies lowest in energy. A third possibility, the 
lambda configuration, occurs when the central state lies highest in energy. 

Unfortunately, I do not have enough space here to present and discuss the 
results of our coherency transfer calculations in stellar atmospheric models for 
all these three-level systems. I have to make a choice: the lambda configuration 
with Ji = 1/2, J2 = 1/2, and J3 = 3/2. My motivation behind this choice is 
two-fold. First, it is precisely the configuration of the Raman scattering pro- 
cess, in which a photon of one frequency is absorbed and a second frequency 
may be re-emitted. Second, this particular three-level atomic system will help 
to understand physically our successful modeling of the "enigmatic" fractional 
linear polarization observed by Stenflo et al. (2000) in the Ca 11 IR triplet (see 
Manso Sainz and Trujillo Bueno 2001a). To this end, the atomic parameters 
of the 1^2 transition have been chosen to be identical to that of the H line 
of Ca 11, while those of the 2^3 transition correspond to that of the 8662 
A line of the Ca 11 IR triplet. The observations were considered "enigmatic" 
especially because this 8662 A line should be, according to Stenflo et al. (2000), 
intrinsically unpolarizable. However, their observations show a polarized peak 
(see also Dittmann et al. 2001). They considered this line "intrinsically unpo- 
larizable" because its upper level (and also the lower level of the H-line) cannot 
harbour any atomic alignment (because J 2 = 1/2 and calcium has no hyperfine 
structure). 

Figure 7 shows the self-consistent results for the unmagnetized case. The 
thin line of the left panel gives the variation with line optical depth of the frac- 
tional atomic polarization of the lower level of the 8662 A line (whose J3 = 3/2). 
The thick line shows the variation of the anisotropy factor at this wavelength (it 
actually shows the variation of —Jq/Jq). We see that depopulation pumping due 
to the 2 <-> 3 transition at 8662 A is able to generate sizeable unequal popula- 
tions among their lower-level sublevels. If one now uses our Eddington-Barbier 
relation (i.e. Eq. 30) and takes into account that for a J/ = 3/2 — > J u = 1/2 
transition W = and i?~0.7, one ends up understanding the computed emergent 
fractional linear polarization shown in the rhs figure. This emergent polariza- 
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tion at fj, = 0.1 is totally due to dichroism in the non-magnetized solar model 
atmosphere (see Trujillo Bueno 1999; page 82). In other words, the generation 
of the observed linear polarization in the Ca n 8662 A line necessarily requires 
a transfer process along the line of sight. This is needed in order to allow the 
Q-component of the absorption matrix (i.e. tjq, which is non-zero because the 
lower level 2 D 3 / 2 is polarized) to make its influence on the emergent linear po- 
larization 0. This is precisely the physical origin of the linear polarization of 
the chromospheric Ca n 8662 A line observed close to the "quiet" solar limb. 

It is important to point out that the calculated line-core fractional polar- 
ization amplitude (Q/J~1.5% as seen in Fig. 7) is an order of magnitude larger 
than the observed one. This "discrepancy" is due to two reasons: (1) because 
the result of Fig. 7 corresponds to a very simple model atom with only three 
levels and (2) because the magnetic field has been assumed to be zero. As shown 
in this volume by Manso Sainz and Trujillo Bueno (2001a) the line-core ampli- 
tude for the zero magnetic field reference case is reduced by about a factor 4 
when a realistic five-level atomic model for Ca n is considered (which includes 
the collisional coupling between the two metastable levels 2 D 3 / 2 and 2 D 5 / 2 ). The 
remaining "discrepancy" is likely due to the Hanle effect of the chromospheric 
magnetic field, as illustrated in Section 7. 

Another interesting example of a group of lines having the same angular 
momentum values and showing a very similar observed behaviour is given by 
the three lines of multiplet no. 2 of Ba n (see the atlas of Gandorfer 2000; 
and note that 82% of the total barium abundance has no hyperfine structure 
splitting). 

6.2. Three lines with a common upper level: the Mg I b lines 

Consider the case of multiplet 3 P° — 3 S, which leads to three spectral lines with 
a common upper level. An interesting example is given by the three lines of 
multiplet no. 2 of Mg I (i.e. the well-known Mg b lines). The Mg 64 line at 5167 
A is a J\ = — > J u = 1 transition (with W = 1 and Z = 0; see Eq. 30), the 
Mg b 2 line at 5173 A is a Ji = 1 -> J u = 1 transition (with W = Z = -0.5), 
and the Mg b\ line at 5184 A is a J\ = 2 — * J u = 1 transition (with W = 0.1 
and Z&0.6). 

Stenflo et al. (2000) considered the Mg b lines "enigmatic" because the 
observed linear polarization amplitudes are similar (see also Stenflo et al. 1983), 
which is in total disagreement with the relative polarization amplitude scaling 
they were expecting from the conventional polarization transfer theory which 
neglects lower-level atomic polarization. 

Which is the true prediction of the conventional polarization transfer theory 
that neglects lower-level atomic polarization? Let us assume that the exact 
fractional atomic alignment of the common upper level at line-core optical depth 
unity along the line of sight is o~q(u) = X > so as to yield, from Eq. (30), 
a positive linear polarization amplitude for the Mg 64 line (i.e. (Q/tybi^X), in 
agreement with the observations (recall that the Mg 64 line has J\ = and that 



Note that the Q component of the emission vector is zero because the upper level cannot 
harbour any atomic alignment. 
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Figure 8. The variation with the corresponding line optical depth of 
the fractional atomic polarization of the lower level of the Mg b\ line, 
of the lower level of the Mg 62 line and of the upper level of the Mg 64 
line. Each of the three optical depths r is measured along the radial 
direction. The solar chromospheric model used is the VAL-C model of 
Vernazza et al. (1981). 



90% of Mg has zero nuclear spin). If we now apply Eq. (30), but neglecting the 
dichroism contribution coming from the atomic polarization of the lower levels 
of the Mg 62 and b\ lines, we find (Q/I)b2~ — X/2 (i.e. a negative polarization!) 
and (Q/I)bi~X/10 (i.e. ten times smaller than that of the 64 line). However, as 
mentioned above, the observations show positive polarization for the three lines, 
and with similar amplitudes. 

What could be the physical origin and explanation of the above-mentioned 
observations? In a previous keynote article (cf. Trujillo Bueno 1999) I suggested 
that one could explain such spectropolarimetric observations if we had (at line- 
core optical depth unity along the line of sight) the following amounts of atomic 
polarization in the metastable lower levels of the Mg 61 and 62 lines: 



• For the lower level of the 62 line — ► <Tq(Z) ~3X. 

• For the lower level of the b\ line — ► ctq(Z) ~ — IX. 
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This "explanation" is based on the direct mapping that Eq. (30) implies: the 
observed polarization amplitudes are giving us directly the atomic polarization 
of the lower and upper levels at optical depth unity along the line of sight ! 

The crucial question now is: are these the atomic polarization values ob- 
tained from self-consistent multilevel scattering polarization calculations in semi- 
empirical models of the solar chromosphere? The answer is affirmative, as indeed 
results from a realistic multilevel coherency transfer calculation for a 19-level 
model atom for Mg I - Mg II (see Fig. 8). The results of this figure show 
that the metastable lower-levels of the Mg I b\ and 6 2 lines are polarized by the 
anisotropic radiation field of typical semi-empirical solar chromospheric models. 
The relative amplitudes of the calculated emergent polarizations of the three Mg 
6-lines at /i = 0.1 agree fairly well with the observations, with dichroism playing 
again a critical role. 

7. The lower-level Hanle effect 

The lower-level 2 D 3 / 2 of the Ca II 8662 A line is metastable. Therefore, from 
Eq. (3), its atomic polarization and the emergent linear polarization have to be 
sensitive to milligauss fields. This is illustrated in Fig. 9, which shows the self- 
consistent values of the /?g-elements for a horizontal field of 10 milligauss^]. The 
assumed horizontal field lies in the xy-plane of the reference system indicated 
in the "north solar pole" example of Figure 1. It is orientated at 45 degrees 
(measured from x towards y). Fig. 9 shows the fractional population imbalances 
(Pq/Po) and the coherences {Pq/pq, with Q = 1 and Q = 2) in percentage. 
Note that the 10 milligauss of the assumed horizontal field are sufficient to 
reduce the population imbalances of the 2 D 3/ / 2 level by a factor 3 approximately 
(compare with the thin solid line of Fig. 7). Moreover, we now have non- 
zero coherences, which are of the same order of magnitude as the population 
imbalances themselves. The emergent fractional linear polarization {Q/I and 
U/I) close to the stellar limb (p = 0.1) and along the x-axis of Fig. 1 (i.e. a 
direction of observation with azimuth x = and p = 0.1) is shown in Figure 10. 
Note that the Q/I amplitude is a factor 3 smaller than the one corresponding 
to the zero magnetic field reference case of Fig. 7, and that U/I is of the same 
order of magnitude as Q/I. 

What is the physical origin of this U/I profile? It is exclusively due to the 
coherences of the metastable lower-level of the Ca II 8662 A line, because t\j = 0, 
while rju^O (see the paragraph after Eq. 25). In fact, it is possible to derive 
Eddington-Barbier relations for the Q/I and U/I emergent polarizations due 
to the Hanle effect. This can be done along similar lines to those which lead to 
Eq. (30) for the unmagnetized reference case. Particularizing the Eddington- 
Barbier relation for U/I to the direction of observation chosen to obtain the 
exact results of Fig. 10 (i.e. x = 0) one obtains: 

U/I « Z y/l[jl- /u 2 Im(a 2 ) + plm(4)], (37) 



These density-matrix values were calculated numerically for the same three-level model atom 
mentioned in Section 6.1. 
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Figure 9. The /jg-elements for a horizontal field of 10 milligauss 
(solid lines) and of 10 gauss (dashed lines). Note that Re^l/Po = 0. 




Figure 10. The emergent Q/I and U/I of the 8662 A line corre- 
sponding to the two cases of Fig. 9. The simulated observation is 
along the x-axis of Fig. 1 (i.e. at % = 0, but at /j, = 0.1) and for an 
orientation of the assumed horizontal magnetic field as indicated in the 
text. 
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where o\ = Pi(0/Po(0 an< ^ °2 = pKO/PoCO are the fractional coherences of the 
lower level of the 8662 A line at the spatial point along the line of sight where 
the optical depth at the selected frequency is unity. Note that i?~0.7, with Z 
the symbol introduced after Eq. (30). Therefore, a positive detection of a U/I 
signal close to the limb (fj, = 0.1) would be a clean observational demonstration 
of the existence of metastable-level coherences. 

Let us now consider what happens when the intensity of the same horizon- 
tal field is 10 gauss. As indicated by the dashed-lines of Figs. 9 and 10, the 
/?g-elements and the emergent linear polarization are modified. The population 

imbalances of the 2 D 3 / 2 level are not destroyed and some of its coherences re- 
main.^. Note also that the emergent Q/I in this saturated gauss regime still has 
a sizeable amplitude (Q//~0.4%), which is due to the atomic polarization of the 
metastable 2 D 3/ / 2 lower-levef^l. In other words, the atomic polarization of long- 
lived atomic levels in the solar chromosphere may still be sufficiently significant 
(even after the partial destruction caused by a 10 gauss purely horizontal field 
with a deterministic or random azimuth!) so as to be able to lead to emergent 
Q/I signals in the observable range. In a forthcoming publication we will show 
in detail that a consistent explanation of the observed Q/I in the three lines 
of the Ca n IR triplet does not necessarily requires to constrain the magnetic 
field to be either extremely low (with B<W milligauss) or vertically orientated 
(but with an intensity in the gauss regime). Note that I am not saying that the 
chromospheric magnetic field is predominantly horizontal. I am basically point- 
ing out that some of the "enigmatic" features of the second solar spectrum are 
due to metastable-level atomic polarization and that the amount of this atomic 
polarization in topologically complex magnetic field scenarios with intensities in 
the gauss regime can still be sufficiently significant so as to lead to the observed 
Q/I polarization amplitudes. 

Finally, it may be found useful to distinguish between five Hanle effect 
regimes in relation with the polarization of a spectral line whose lower level 
is the ground or a metastable level. These Hanle regimes can be established 
according to the values of T[ and T u and by noting that, for a given magnetic 
field intensity, Ti 3> T u (because the lower-level lifetime is about two orders of 
magnitude larger than the upper- level lifetime). As mentioned earlier, in this 
section we are using the reference system whose z-axis is directed along the 
radial direction of the star:[^] 

• Regime 1: T t -C 1 and r u < 1 

In this regime the atomic alignment of the lower level of the line transition 
of interest is typically of the same order of magnitude as the upper-level 



For the case of a microturbulent and isotropic magnetic field or for a microstructured horizontal 
field with random azimuth, only the population imbalances would remain (see the dashed- 
dotted line of Fig. 3 and the Appendix of the paper by Trujillo Bueno & Manso Sainz 1999). 

We still find Q/I w 0.1% when we use, instead of a three-level model, a realistic 5-level atomic 
model including the collisional coupling between the metastable levels 2 D 3 / 2 and 2 D 5 / 2 . 

See Landi Degl'Innocenti (1985) and Section 5.2 for a discussion of Hanle-effect regimes from 
the point of view of the magnetic field reference frame. 
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alignment (i.e. | cr§ (Z) | ~ |<7q(u)|). We find here the maximum possible Q/I 
amplitude. 

• Regime 2: « 1 and T u <^ 1 

In this regime the emergent polarization is sensitive to milligauss fields. 

• Regime 3: T t > 1 and r u < 1 

In this regime we have |crg(Z)| <C |ctq(u)|. 

• Regime 4: F[ » 1 and T M ~ 1 

In this regime the emergent polarization is sensitive to magnetic fields in 
the gauss regime. 

• Regime 5: > 1 and T M > 1 

We again find that |crg(Z) | ~ |<7Q(it)|, but with lower individual values for 
|oq(/)| and |o"o(u)| than for the unmagnetized reference case. The emer- 
gent polarization is sensitive only to the field orientation. In this regime, 
one should find relative amplitudes of the linear polarizations observed in 
multiplets (e.g. in the Mg I b lines) similar to those found for the unmag- 
netized reference case. This emphasizes the importance of an accurate 
quantification of the observed polarization amplitudes in order to infer 
correctly the intensity and inclination of the magnetic field. 

8. Concluding remarks 

This review article, besides providing an introduction to optical pumping, atomic 
polarization and the Hanle effect, has advanced some new results which we 
will publish in suitable journals during the following months. Of particular 
interest for the reader who looks at the Sun as a unique physics laboratory is 
the conclusion that the quiet solar chromosphere is a polarized vapor with optical 
properties similar to those of an anisotropic crystal. 

On the other hand, the reader interested mainly in the remote sensing of 
solar and stellar magnetic fields should feci optimistic because the quest for 
understanding the physical origin of the "second solar spectrum" is now becom- 
ing one of the success stories of astrophysics. We had the theory (cf. Landi 
DegPInnocenti 1983). We had the observations (cf. Stenflo and Keller 1997; 
Stenflo et al. 2000). And now we have the self-consistent radiative transfer 
modeling, which is the essential link between theory and observations. In the 
following years rigorous confrontations of observations of weak polarization sig- 
nals with numerical simulations of scattering polarization and of the Hanle and 
Zeeman effects should lead to fundamental new insights in our understanding of 
solar photospheric and chromospheric magnetism. 
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